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Strictly evolutionary approaches to improving the air transport system — a highly com- 
plex network of interacting systems — no longer suffice in the face of demand that is pro- 
jected to double or triple in the near future. Thus evolutionary approaches should be 
augmented with active design methods. The ability to actively design, optimize and con- 
trol a system presupposes the existence of predictive modeling and reasonably well-defined 
functional dependences between the controllable variables of the system and objective 
and constraint functions for optimization. Following recent advances in the studies of the 
effects of network topology structure on dynamics, we investigate the performance of dy- 
namic processes on transport networks as a function of the first nontrivial eigenvalue of 
the network’s Laplacian, which, in turn, is a function of the network’s connectivity and 
modularity. The last two characteristics can be controlled and tuned via optimization. We 
consider design optimization problem formulations. We have developed a flexible simu- 
lation of network topology coupled with flows on the network for use as a platform for 
computational experiments. 


Nomenclature 


Ai The smallest eigenvalue of L 
X 2 The smallest nontrivial eigenvalue of the L 
Xmax The largest eigenvalue of the L 
A Adjacency matrix 

Oij Elements of A; 1 if node i connected to node j; 0 otherwise 

c(G) Clustering coefficient 

D Degree matrix 

di Degree of node i 

G Graph corresponding to a network 

L Network Laplacian 

N The number of nodes in G 

s{G) A measure of the tendency of nodes of similar degree to be connected 
Smax Maximum of s(G) 

Smin Minimum of s(G) 


I. Introduction 

L arge technological systems, such as the Internet and the air transportation system, have never been 
actively designed in their entirety; instead, they evolved gradually, in response to demands. Accommo- 
dating the projected near-future threefold increase in demand for air transportation, as well as the growing 
complexity of airspace, via strictly evolutionary changes is problematic. There appears to be a clear need 
for active and rigorous design methods. 
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Active rigorous design necessitates identification of design variables, as well as objectives and constraints 
whose functional dependence on the independent variables has been established analytically and experimen- 
tally. In other words, active design requires predictive models. Although design of systems such as air 
transportation may never be as amenable to active design as are complex physical artifacts (e.g., airplanes, 
engines), we conjecture that a significant portion of active design is possible with appropriately developed 
variables and metrics. 

Air transportation is a highly complex network. The past several years have seen a surge in the study 
of both the structure (see, e.g., Newman^ for an overview) of complex networks and the dynamics (e.g., 
Strogatz") of flows on complex networks. The effect of network structure on the dynamics is the starting point 
in our investigation aimed at deriving quantifiable functional relations between global transport metrics (e.g., 
throughput, delays, capacity) and locally controllable structure (e.g., connectivity) of transport networks. 
Of particular interest is computational controllability of the local variables expressed as an optimization 
problem. 

Recent groundbreaking work in a variety of domains (statistical physics, biology, Internet router design) 
derives models that link static topological network structure variables, such as degree distribution, clustering 
coefficient, and average shortest path distance, with the dynamic performance of flows over networks (e.g., 
synchronicity, throughput, congestion). In particular, Atay et al.,’’"^ Donetti et al.,^ Motter et al.®’^ have 
begun to establish explicit functional dependences between topology and system dynamics on networks. 
Other measures, such as assortativity and s(G) (the tendency to connect with nodes of like degree) has been 
explored by other authors, e.g., Li et al.®’® 

While there is a superficial similarity between the systems under study in other domains and the air 
transport system, there are also substantial differences. For instance, the variation in the Internet packet 
size is not nearly as great as the variation in the plane capacity; the connections (routes) among nodes of the 
airport network are not as hard wired as for the Internet. Moreover, the assumptions made in the physics 
literature are not applicable to the transport network. For example, Atay et al. assume that all nodes 
are identical and conform to a generic discrete time equation. These assumptions are needed to simplify 
the problem in the initial stages of research, but a question remains whether the derived models would be 
sufficiently applicable to air transport. 

Having conducted initial computational experiments, we believe that there is much to be learned from the 
research efforts in other domains and, although the models will not carry over one-to-one to air transport, 
the models from other domains can be developed further to accommodate transport assumptions and needs. 
In this paper, we summarize our results to date: a computational study of network topology metrics; 
construction of a test platform; and formulation of optimization problems. We will call our problem air 
transport network (ATN) optimization. 

II. Spectral Graph Properties and Dynamic Processes 

Several authors®’ ^ have recently examined network dynamics as a function of network topology and 
showed that different constrained topology optimization problems, such as maximizing synchronization or 
node proximity, lead to optimal topologies that, although not identical, share common features. Donetti et 
al.® call these optimal networks entangled. The result appears relevant to optimization of transport networks. 
In this section, we review some of the basic concepts. 

Let G be a graph (network) with N nodes. Associated with G is the adjacency matrix A = [a^], where 
Oij = I if nodes i and j are connected and 0 otherwise. The degree matrix D = [di\ is a diagonal matrix 
of node degrees. The Laplacian matrix L = D — A fully captures the connectivity structure of G. The 
eigenvalues of L satisfy 

‘^dmax ^ ^ ^ ^ — 0 ; 

where d^ax is the maximum degree of G, and the eigenvector (!,...,!) corresponds to the trivial eigenvalue 
The first nontrivial eigenvalue A2 is known as the spectral gap or algebraic connectivity. 

Donetti et al.® give an illustrative example that foreshadows the topological significance of the spectral 
gap for networks dynamics. If a network consists of a number of disconnected subgraphs, its Laplacian is 
block-diagonal and the multiplicity of the trivial eigenvalue equals the number of disconnected subgraphs. 
Connecting the subgraphs weakly introduces small eigenvalues with nearly constant corresponding eigenvec- 
tors. This feature (small spectral gaps) provides criteria for graph partitioning in well-known algorithms. 
Intuitively, small A2 values imply the existence of well-defined modules that can be disconnected by cutting 
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a small number of links, while large A 2 values point to unstructured (entangled) graphs. 

Several authors studied synchronizability of diffusive processes on networks, considering a general 

dynamical process 

N 

Xi = F{xi) - a'^LijH{xj),i = ( 1 ) 

i=i 

where Xi are dynamical variables, F is an evolution function, is a coupling function, and cr is a coupling 
constant. Although diffusive processes are known to have synchronous states, the question is, under what 
conditions these states are stable. A linear stability analysis reveals that synchronized states are more 
stable for smaller Since the variability of the maximum eigenvalue is bounded for graphs under 

consideration, increasing stability of synchronized states amounts to maximizing the spectral gap A 2 . Other 
authors ’’ have used the spectral gap as an indicator of synchronization for discrete processes as well. 

The normalized Laplacian, L' = D~^L, and its eigenvalues {A'} also play an important role, especially 
in the study of random walks, a subject relevant to propagation of traffic through networks. Large spectral 
gaps increase the rate at which random walks move and disseminate. A class of graphs with large spectral 
gaps, known as Ramanujan graphs, is described by Donetti et al. ’ These graphs are regular, have a vanishing 
clustering coefficient, a small average shortest path distance and a large girth. Since air transport networks 
do not have this structure, general model (1) is not applicable. However, changing the coupling constant 
to and thus normalizing the effect of the neighboring nodes (in turn, increasing the relevance to traffic 
networks), results in an optimal topology when the normalized spectral gap is maximized. 

In studying networks metrics, we also examine the function, s(G), developed by Li et al.®’*^ This metric 
distinguishes among network structures of a fixed degree distribution, based on the tendency of nodes of 
similar (or dissimilar) degrees to be connected. It serves as an indicator of performance in, e.g., router 
design. The metric is defined as 

s(G)= ^ d,xdj, (2) 

{i,j)^Links 

where di is the degree of node i. The metric is often normalized as 


S{G) 


<G) - 


<G) - 


where Smax and Smin are the maximum and minimum values, respectively, of the sum taken over all possible 
connected graphs for a fixed degree distribution. Since it turns out that Smin is difficult to compute in 
practice, the metric is often approximated by 


S{G) 


s{G) 

s(G) Smax 


Li et al.® term graph realizations with large values of s(G) scale-free and graph realizations with small 
values of s(G) scale-rich. Consequently, in scale-free graphs high-degree nodes are more likely to be adjacent 
to other high-degree nodes, while in scale-rich graphs high-degree nodes are more likely to be adjacent to 
low-degree nodes. 


III. Transport Network Design Formnlations 

Entangled topologies optimal for such objectives as synchronization and maximum proximity of nodes 
bear some resemblance to topologies of transport networks derived by Guimera et al. for transport with 
congested hubs. Intuitively, this makes sense: if hubs are congested, then the best way to disperse the load is 
by providing a massively connected network that is difficult to separate by cutting a few links. This supports 
the use of the spectral gap as one of the objectives for ATN design. 

Our initial experiments with the spectral gap^^ further support the notion of its usefulness as an indicator 
of performance. However, many questions remain open. As we mentioned, the assumptions underlying the 
construction of optimal entangled topologies are not, in general, applicable to realistic transport networks. 
For instance, transport networks are not regular; their degree distributions are not necessarily fixed. We 
need to introduce assumptions that will eventually allow us to model more practical conditions. 
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As a start, consider that the spectral gap of the entire network is too global a measure. We do not have 
the freedom of rewiring the network completely, especially considering the many players in the network (e.g., 
the air carriers). Moreover, even in the presence of congestion, doing away with hubs completely may not be 
an option because hub-like structures have a number of important benefits, such as the small world property 
(minimal hops) and robustness with respect to random disturbances. The optimal answer likely lies in a 
mix of modularity and entanglement. Is there a way to pose the graph optimization problem that preserves 
some degree of global modularity, while allowing for local entanglement within modules? 

Let Go, Gi, . . . , G„ be a decomposition of G, as follows. Go represents the subset of the network shared 
by all other subsets. This subset must maintain some degree of modularity. Let Vi be the set of nodes 
(vertices) of graph i. Gi, . . . , G„ are subgraphs allowed to entangle locally. All subsets must maintain a 
fixed s(Gi) and degree distribution 

«w = v S 1. 

* VjGVi\dj<d 

where Ni is the number of nodes in subgraph i. 

We will use problem formulations analogous to those used for multidisciplinary design optimization 
(MDO).^* However, due to a different nature of ATN (discrete and, eventually, mixed), the analogy to MDO 
formulations is purely structural. For now, the number and location of nodes will be fixed (as in the present 
airport network). Consider an example from a range of formulations being studied. 

In this formulation, we employ bilevel optimization with no feedback. Let Q denote simple connected 
graphs. The outer problem will maximize the spectral gap of the global graph Gq, i.e., we solve approximately 

minimizeGg X2(Gq) 
subject to s(Go) = G^“ 

Po{d) = 

Go&G, 

G G 

where, Gi and G 2 are constants. Then the subproblems will minimize their spectral gaps to increase 
connectivity, subject to maintaining the inter-subgraph connections fixed. That is, for each subproblem, we 
solve approximately 

maximizeci A2(G^) 

subject to s(Gi) = Gf ‘ , , 

P,{d) = 

Gq G Q remains fixed. 

Formulation (3-4) should maintain connections among the modules, while maximizing the fluency of the 
flow inside the modules. Other authors^^ have considered other notions of partitioning the network into 
local and global subsets. 


IV. Test Platform 

Many sophisticated airspace simulation models have been developed in recent years. They include 
trajectory-based models, such as ACES^'’ and FACET;^' models based on network flow with physical analo- 
gies;^^ empirical models;^^ models based on Bayesian networks."^ See the last reference for a more complete 
list of models. The existing models operate under a variety of assumptions and predict various effects, such 
as congestion and delays. Despite the availability of a number of models, we have decided to develop our own 
research test platform. We do this for a number of reasons. To derive functional dependences in the system, 
we need a completely flexible definition of the system components. This implies control over the source code. 
Computational expense of some of the existing detailed airspace models is another consideration. Finally, 
we can use data generated by the existing models as a load on our system. 

The Airport Network Simulation Program (ANSP) was designed as a first step in the development of 
an optimization framework able first to investigate and, we hope, eventually to model and design complex 
transport networks in the context of the total transportation system.^® Version 1 of ANSP is a discrete-event 
simulation, with events separated by discrete time intervals and processed chronologically. A global event list 
drives the simulation by coordinating interactions between the simulation’s two primary structures: nodes 
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(airports) and planes. While events are currently limited to arrivals, departures and vector updates, ANSP’s 
flexibility allows new events to be added easily. The target is an adaptive, gradually expanding simulation 
that will assist in numerical experimentation with transportation system concepts and design methods. 

At completion, each simulation computes a number of metrics, currently including the total holding time 
for each flight, the average number of planes in queue for each airport; total number of hops, time in air, 
and distance traveled for each flight. 

V. Computational Experiments 

To-date, we have experimented extensively to observe the relations among various network metrics, such 
as the spectral gap, assortativity (as s(G)), cluster coefflcient, graph eccentricity and diameter, and others. 
The experiments confirmed that the spectral gap should prove useful, at some level of detail, for ATN design. 
They also raised questions about the correlation among the metrics. To answer these questions, we must 
transition from the purely static structure experiments to experiments that include trafflc flow dynamics. 

In this section, we give an example of static experiments with the spectral gap and report on the initial 
test of the ANSP platform that will allow us to introduce dynamics on the network. ANSP tests were conducted 
using actual 1990 U.S. Air Transport data.^^ 

V.A. Experiments with the spectral gap 

We experimented with two types of abstract, 100-node, simple, undirected networks (no loops or multiple 
links): preferential attachment and geometric. Preferential attachment graphs are generated following the 
approach given by Barabasi,^' where a network is grown by adding nodes and edges; for each node added, 
m edges are added preferentially, based on the current degree distribution. We chose m = 2 to facilitate 
visualization. 

Geometric graphs are generated by randomly selecting 100 points (r, 9) with values of r G [0, 1] and values 
of 0 G (0, 360]. Edges exist between pairs of points if the Euclidean distance is less than a specified threshold 
(in our experiments thresholds between 0.17 and 0.25 were used). If the resulting graph is connected, it is 
kept; otherwise it is rejected and the process begins again. 

The graphs plotted in Figures 1 and 2 were constructed by first generating a random instance of the 
particular graph class — geometric in Figure 1 or preferential attachment in Figure 2. Next, a simple tabu 
search^® heuristic was called to minimize or maximize the spectral gap, while keeping the degree distribution 
and s(G) fixed. Allowable re-wirings are pair-wise edge interchanges that preserve the degree distribution 
and s(G). Briefly, the tabu search checks to see if the move is acceptable, that is, if the move is improving 
and not tabu, or improving and tabu but leads to the best observed value of A 2 (aspiration criterion). 

Figures 1 and 2 display networks with respect to the reciprocal of the eccentricity of each node u. The 
eccentricity of u is its maximum (shortest path) distance. The graphs are generated by socnetv^. The goal 
of the plots is to uncover any qualitative differences between the graphs with small and large values of the 
spectral gap. Nodes with equal eccentricity values are plotted on the same (dashed line) circles. The circles 
with larger radii have larger eccentricity. Consequently, nodes near the center have shorter longest paths. 

Qualitatively, when A 2 is small, the patterns are less organized, the eccentricity plots in Figures la and 
2a are more dispersed and consist of many rings of constant eccentricity. The eccentricity plots with larger 
A 2 are more organized, with few rings of constant eccentricity. We refer the reader to an earlier paper for 
detailed description of the experiments and other metrics. 

V.B. Experiments with U.S. 1990 air transport data 

The functionality of ANSP was tested by comparing simulation statistics for the original network and several 
rewirings produced by tabu search.^* It is important to note that at this stage, ANSP alone cannot capture 
the true system complexity and dynamics; these computations are used purely to verify the performance of 
the platform, rather than to draw any definitive conclusions about the functional dependences. 

U.S. 1990 airport activity statistics^® were used to generate an airport network with 336 nodes (airports). 
The network was then optimized, using four runs of tabu search. The four resulting networks were “optimal” 

®'The source code and documentation can be found at http://socnetv.sourceforge.net/. 
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(a) A 2 = 0.009, c{G) = 0.426 (b) A 2 = 0.314, c{G) = 0.297 

Figure 1. Geometric graphs: 100 nodes, s(G) = 0.971, fixed degree distribution 



(a) A 2 = 0.006, c{G) = 0.210 (b) A 2 = 0.365, c{G) = 0.101 

Figure 2. Preferential attachment: 100 nodes, s{G) — 0.716, fixed degree distribution 
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Figure 3. 


Average holding time and distance traveled for 1990 U.S. airport network vs. its rewired counterparts 


with respect to one of the four objectives (smin, Smax, ^2min, and \2max), subject to preserving the node 
degree sequence. Again, since tabu search is a heuristic, true optimality is not guaranteed. 

We now had five networks: the original one and four networks optimal with respect to one of the four 
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Figure 4. Average number of hops for 1990 U.S. airport network vs. its rewired counterparts 


metrics. The original airport network, s{g) « 104, 276, 109 and A 2 « 0.339. The tabu search algorithms 
produced four networks with the following heuristically optimal metrics: 

Smin ~ 90, 988, 359 Smaa; ~ 109, 205, 715 A2mm ~ 0.132 A2ma£c ~ 0.382. 

Now the five networks were imported into ANSP, which ran 100 independent trials, each generating 400 
random origin-destination routes. We chose arbitrary large numbers of trials and routes to compute statistical 
averages of the total holding time (time in a holding pattern), total distance traveled, total time in air, time 
average number of planes in queue (arrival and departure) at each airport, and total hops for each flight 
were computed for each trial. The results were compared across all networks and are shown in the plots of 
Figures 3-5. (Time is in hours.) 

In this experiment, we assign superior performance to lower numbers for holding time, distance traveled, 
time in air, and hops. A qualitative analysis shows that the original network and the \ 2 max rewiring exhibited 
the best overall performance. Both networks have relatively large A 2 values, and hence high connectivity, 
relative to other networks with the same node degree sequence. The \ 2 min network was a low performer. 
We conjecture that low connectivity and the reliance on indirect shortest paths for flight routes may be a 
contributing factor. The Smax network exhibits similar connectivity. However, the low performance of the 
Smin network cannot be attributed directly to its spectral properties, as its spectral gap is large. Further 
analysis shows that the greatest performance detractor for the Smin network is the reliance on major hubs. 
This may appear counterintuitive, since a smaller s(G) typically corresponds to a less hub- like topology. 
However, the preponderance of east-west flights in the network appears to magnify the Smin system’s use of 
hubs, specifically the Chicago O’Hare International Airport (ORD). After 100 simulations, approximately 58 
flights per trial passed through ORD in the Smax network and 95 for the Smm network. ORD was consistently 
the most loaded airport in each network, which is a clear contributing factor to large average holding time 
and average time in air for the Smin network. A tentative conclusion is that, while the spectral gap and 
s{G) are informative for linking topological and dynamical properties, further assumptions and modeling 
elaborations are needed to account for realistic dynamic behavior. 

A second test was conducted on a smaller network of 318 nodes (airports) that did not include Alaskan 
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Figure 5. Average time in air and in queue for 1990 U.S. airport network vs. its rewired counterparts 


airports. Some of the metrics are summarized in Tables 1-3. Table 1 contains metrics from the original 
network, Table 2 from the network rewired for Smax and Table 3 from the network rewired for Smin- Note 
that the removal of Alaska increased the value of the spectral gap, compared to that of the 336 node network. 
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"^ 2 , min 

A2,max 

"^ 2 .original 

A 2 

0.168907 

0.583535 

0.499723 

Xn 

199.020485 

199.020809 

199.020596 

s(G) 

92613477 

92613477 

92613477 

c(G) 

0.6160 

0.6145 

0.6164 


Table 1. 1990 Regular 



'^ 2 , min 

A2,max 

^ 2 , original 

A 2 

0.172926 

0.580941 

0.507063 

Xn 

199.011316 

199.011787 

199.011806 

s{G) 

95354201 

95354201 

95354201 

c(G) 

0.6575 

0.6575 

0.6575 


Table 2» 1990 



"^ 2 , min 

A2,max 

'^ 2 , original 

A 2 

0.970814 

0.971974 

0.971191 

Xn 

200.074690 

200.084046 

200.079116 

s{G) 

80767804 

80767804 

80767804 

c(G) 

0.4186 

0.4185 

0.4185 


Table 3. 


1990 s 


min 


This stands to reason as removing a “module” lead to a more entangled network. 

VI. Concluding Remarks 

The present state of air transportation networks necessitates developing active and rigorous design 
methodologies to augment the traditional evolutionary approaches. We investigate the applicability of net- 
work models originating in statistical physics, Internet router design, and biology to air transport design. 
Although some of the modeling assumptions do not hold, the initial results are sufficiently interesting to 
warrant an effort in adapting the models for use in air transport analysis and optimization. We are cau- 
tiously optimistic that introducing increasingly realistic formulations will allow us to progress from abstract 
networks to models that approximate transport networks to a greater degree of fidelity. 
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